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We calculate the backflow current around a fixed impurity 
in a Fermi liquid. The leading contribution at long distances 
is radial and proportional to 1/r 2 . It is caused by the current 
induced density modulation first discussed by Landauer. The 
familiar 1/r 3 dipolar backflow obtained in linear response by 
Pines and Nozieres is only the next to leading term, whose 
strength is calculated here to all orders in the scattering. In 
the charged case the condition of perfect screening gives rise to 
a novel sum rule for the phase shifts. Similar to the behavior 
in a classical viscous liquid, the friction force is due only to the 
leading contribution in the backflow while the dipolar term 
does not contribute. 
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The calculation of backflow in liquids is one of the stan- 
dard problems in hydrodynamics, determining e.g. the 
Stokes friction in a classical viscous liquid or the prop- 
erties of rotons and impurities in superfluid Helium 4 [^J . 
For the case of a Fermi liquid, the backflow around a 
slowly moving massive impurity is discussed in the clas- 
sic text by Pines and Nozieres || . Within linear response 
they show that the backflow is proportional to the density 
response function and dipolar in character. As pointed 
out by these authors, the dipolar form may be derived 
from a rather simple geometrical argument: Indeed the 
backflow current outside the impurity should have zero 
divergence being a stationary flow and zero curl because 
the perturbation is longitudinal. The only vector func- 
tion obeying both conditions, however, is a dipole. For 
a neutral Fermi liquid the strength of the dipole is given 
by the compressibility times the Fourier transform of the 
interaction at zero momentum. In the charged case the 
dipolar backflow has a universal amplitude. This is a re- 
sult of perfect screening which requires that the backflow 
identically cancels the longitudinal part of the impurity 
current |^]. 

In this letter we reconsider the backflow problem in a 
Fermi liquid, going beyond the linear response treatment. 
Starting with the simple case of a nonintcracting Fermi 
gas, we show that the leading term at long distances is 
not the dipolar backflow but a radial contribution decay- 
ing like \jr d ~ x in d dimensions (d = 2, 3 in the follow- 
ing). It is proportional to the impurities transport cross 
section and thus is not contained in a linear response 
calculation where the interaction only appears linearly. 
The novel term has nonzero curl and is directly related 
to the asymmetry in density around localized scatterers 
in the presence of a finite transport current, discussed 



long ago by Landauer Q. We also calculate the famil- 
iar dipolar backflow to all orders in the scattering po- 
tential and discuss the modifications of our results for 
interacting Fermi liquids. It is found that the condition 
of perfect screening entails a sum rule for the scattering 
phase shifts which is similar to, but different from the 
one by Friedel ||. Finally we determine the systematic 
force exerted on the impurity by the moving fcrmions. It 
is shown that only the leading l/r d_1 term of the back- 
flow current contributes to the force, a situation which is 
completely analogous to that in a classical viscous liquid. 

Let us consider a fixed scattering center at the origin 
which is characterized by a spherically symmetric inter- 
action potential V{x). In the frame where the impurity 
is at rest, the Fermi system is flowing past with asymp- 
totic velocity v ^ 0. The unperturbed current density 
is therefore j(x)\o = nv, with n the equilibrium number 
density. Due to scattering off the impurity, the actual 
current density j(x) differs from nv by a backflow cur- 
rent Sj(x). To lowest order in v the Fourier transform 
Sj(q) of the backflow is of the form 



Sj(q) = Hq) [(qv)q - v\ 



(i) 



where q is the unit vector in the direction of q. Indeed the 
vector in Eq.(|lj) is uniquely determined by the require- 
ment that it is linear in v and the zero divergence condi- 
tion q-Sj(q) = due to the stationarity of the flow. For 
small velocities the backflow pattern is thus completely 
determined by the scalar function h(q). As pointed out 
above, a treatment of the interaction potential V(x) in 
linear response gives rise to a dipolar backflow which is 
characterized by lim^o h(q) = ho. The associated di- 
mensionless constant ho is equal to dn/dfi ■ V(q — 0) in 
the case of a neutral Fermi liquid || . Here the compress- 
ibility dn/dfx is just the q — > limit of the general density 
response function x(<z)- F° r an impurity with charge Z, 
the product x(q) ■ V(q) is replaced by Z [e' 1 (q) — l] with 
e(q) the static dielectric constant ||. As a result of the 
perfect screening condition e^ 1 (q — > 0) = 0, this leads to 
a universal value h$ — — Z for the strength of the dipolar 
backflow in the charged case. 

In order to discuss the generalization of these results 
beyond linear response, still keeping the asymptotic ve- 
locity v small however, we start by considering a nonin- 
teracting Fermi gas. In this case the backflow can be cal- 
culated analytically from the single particle eigenstates 
ipk{x), which are the exact outgoing scattering states in 
V{x). Indeed describing the finite asymptotic current nv 
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by a shifted Fermi distribution /(eg_ m ^ ?i ) for the in- 
coming momenta fc, the total fermionic current density 
at zero temperature and to linear order in v is given by 



i(z) = 



„d-l 



(2tt)° 



dQ. k k-vlm[ipl{x)V x ^ k {x)\\k=k F ■ (2) 



Here dfl k denotes an integration over the directions of the 
unit vector fc, while the magnitude k = \k\ is fixed at the 
Fermi wave vector kF- Thus at T = and to linear order 
in v the backflow is completely determined by the exact 
scattering states right at the Fermi surface. Clearly the 
behavior of j(x) at arbitrary distances depends on the de- 
tails of tpk{x)- For large distances, however, it is sufficient 
to know the asymptotic form of the scattering states. In 
order to obtain the first two leading contributions to j(x) 
it is necessary to expand 
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(3) 



to order 1/r 2 in three dimensions. Here / is the standard 
scattering amplitude while the coefficient of the 1/r 2 con- 
tribution is given by 

oo 

h = ^2 X)( 2/ + l W + ^ sin<5 ' P ^ ■ *) (4) 

2=0 

with phase shifts di and the usual Legendre polynomials 
Pi. This result is obtained by a straightforward asymp- 
totic expansion of the free particle solutions with given 
angular momentum I. In two dimensions the correspond- 
ing result turns out to be 
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(5) 



with amplitudes / and ji which are not given here ex- 
plicitely. It is now straightforward to insert the asymp- 
totic behavior of the scattering states into our expression 
(2) for the current. Apart from the trivial term k which 
accounts for the background current density nv, the lead- 
ing contributions to Im[-0* V^] \k=k F obviously arise from 
the square kpx\f\ 2 /r d ^ 1 of the outgoing wave and the 
two interference terms linear in /. Now expi(fcr — kx) is 
asymptotically proportional to a (5-function S(Tlk — Q x ) 
which singles out the forward direction k = x. Using the 
optical theorem it is straightforward to show that the 
leading term to the backflow is given by 



Sj(x) 
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T (xv)x +0(r- d ) 



(G) 



with <j tr = J dflk(l — kx)\f\ 2 the standard transport cross 
section. Obviously the contribution (6) is a purely radial 
current which vanishes in the direction perpendicular to 
v (see Fig.l). It has vanishing divergence as it should, 
but finite curl. 
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FIG. 1. a) Radial backflow current —(xv)x/r in d — 2. 
The impurity sits at the center with incoming current from 
the left, b) Dipolar backflow (v — 2(xv)x) /r 2 in d = 2 for a 
repulsive impurity at the center. The direction of the flow is 
reversed in the attractive case, in contrast to a). 

In order to understand its physical origin we consider 
the current induced part 5n{x) of the density modula- 
tion which is caused by the scattering off the impurity. 
As predicted by Landauer |Q this modulation asymp- 
totically has the form 8n(x) ~ — (xv)/^^ 1 of a dipole 
potential. Comparing the exact expression obtained for 
5n(x) in a scattering theory calculation with our re- 
sult (6), it turns out that at T = and to linear order in 
v the asymptotic backflow current is simply given by 



5j(x) = vpSn(x) ■ x 



(J) 



The leading term in the backflow is thus directly propor- 
tional to the current induced density change Sn(x) which 
is positive in front and negative behind the scatterer, in 
agreement with the intuitive picture developed by Lan- 
dauer Q. As a result, the sign of this contribution to the 
backflow remains unchanged upon going from a repulsive 
to an attractive potential V(x). This is in contrast to the 
dipolar contribution 



Sj{x)\di P 



ho d(xv)x 



2n{d-l) 



(8) 



which is only the next to leading term in an asymptotic 
expansion of j(x). Using (3) and (5) a straightforward 
but rather tedious calculation indeed gives a contribution 
to Sj(x) of the form (8) with strength 
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h = — - ReVo V"c( sin (5 ; sin(5 ;+ i sin (5i - S i+ i) 

on 7r f— ' 

r (=0 



(9) 



Here q = 2/ + 1 or (I + l) 2 in d = 2 or d = 3 respectively, 
while 



y = / f«- kI y(f)4(f) 



(10) 



is essentially the exact forward scattering amplitude. 
Since V reduces to V(q = 0) in Born approximation, 
the first term in (9) is the obvious generalization of the 
linear response result of Pines and Nozieres to arbitrary 
order in the scattering potential. It is convenient to ex- 
press also this contribution in terms of the phase shifts 
Si via 



h(P = ~-^^ a i sin 2(5; 

1=0 
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(11) 



where a; = 2 — S^a or 2/ + 1 in d = 2,3. In addition to h ^ 
there is a second contribution which is at least of order 
V 3 . It arises from the interference term kpx/r d -2Re f*fe 
between the first and next to leading contribution to the 
outgoing wave. The additional term is odd under <5; — * 
—Si as is the first one, but vanishes in the case of s-wave 
scattering only. 

Including both the radial and dipolar contributions to 
the backflow current, the scalar function h(q) defined in 
(1) has the general form 



lim h(q) = 

q ^o w q 



ho + ... (12) 
with h-i = (2, 37r/4)n<7 tr in d — (2,3). Since 



0tr 



k d-l 



h sin 2 (5i - 5i+i) 



(13) 



with bi = 1 or 7r(/+l) in d = 2,3, both leading coefficients 
and ho can be expressed completely in terms of the 
density n ~ k F and the scattering phase shifts Si. In 
the particular case of d = 3 and s-wave scattering with 
scattering length a, we have = (kpa) 2 kp/2 and ho = 
kpa/n. With kp as the typical scale for q, this shows 
that the strength of the leading radial term is then a 
factor kpa <C 1 smaller than the dipolar contribution. 
Nevertheless at long distances it is always the radial term 
which dominates. 

In a Fermi liquid the interacting state develops adi- 
abatically from the noninteracting one. The resulting 
quasiparticle states in a local and centrally symmetric 
potential are characterized by phase shifts Si . In general 
these are functionals of both the energy and the quasipar- 
ticle distribution rik(x) B. By Galileian invariance the 



asymptotic distribution is again a Fermi sphere shifted by 
5k = mv/h, with m the bare mass. At T = the energy 
is fixed at ep and there are no collisions other than with 
the impurity. Moreover since the deviation of from 
equilibrium is already linear in v, we may neglect the de- 
pendence of 5i both on energy and on nu- The resulting 
values of Si then define an effective force Fk on the quasi- 
particles which appears in the corresponding transport 
equation ||. In a fully quantum mechanical treatment 
of the Wigner function rik{x) the associated local parti- 
cle current must then be equal to jk{x) = ^Imip k 'V x 4>k 
where tpu are the exact scattering states in an effective 
potential with phase shifts Si. As was shown above, h-\ 
and ho can be expressed completely in terms of kp and 
the scattering phase shifts Si. The generalization of our 
results to the interacting case is therefore rather obvious. 
Indeed since kp is unchanged one only needs to replace 
the phase shifts by those for quasiparticles. The general 
form of the backflow as determined by (1) and (12) thus 
applies also in the interacting case, however with renor- 
malized parameters and h . For a charged impurity 
in an electron liquid, the perfect screening condition must 
hold to all orders in V . As we have seen, this implies a 
universal dipolar backflow characterized by h^ = — Z for 
an impurity with charge Z. Since ho is completely deter- 
mined by the Si via (9) and (11), perfect screening gives 
a nontrivial condition on the scattering phase shifts at a 
charged impurity. In the limit <5; <C 1 it reduces to the 
well known Friedel sum rule Q which fixes the number 
of bound states. The novel sum rule shows that even 
for Z = 1 no purely s-wave scattering potential can ac- 
count for the backflow in the charged case. Regarding the 
dominant radial contribution, the transport cross section 
appearing in the coefficient has to be replaced by its 
value for the screened potential V(q)/e(q). In contrast 
to ho the strength of the radial backflow is therefore not 
universal. 

Finally we calculate the systematic force F due to the 
transfer of momentum between liquid and scatterer. In 
the context of electromigration theory this is known as 
the wind force M. Taking the gradient of the interac- 
tion energy with respect to the impurity position, it is 
straightforward to see that 



F = — J d d xn{x)V x V . 



(14) 



Clearly at zero current v = this force vanishes al- 
though the fermion density is not uniform even in this 
case. Therefore only the current induced density change 
Sn(x) contributes to F. For simplicity we consider again 
a Fermi gas at T = with scattering states \k+ >. To 
lowest order in v the current induced density then has 
the asymptotic behavior (d = 2,3) 



8n{x) = — 



1 
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(15) 
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which shows that the two leading terms in vpSn(x) are 
identical with the radial component x ■ Sj(x) of the back- 
flow. The associated total force can be written as 

/c d_1 f » m 
F = -Z-rj J dn k k-vj-<k+\- V x V\k+ > \ k = kp 

(16) 

similar to (2). Now the relevant matrix element of 
V X V between the exact scattering states is equal to 
2eF&tr(kF) ■ k. Thus (16) immediately gives a conven- 
tional friction force F = —7]fv with r\p — hkpntjtr @- 
The fermionic friction coefficient t\f is proportional to 
the transport cross section which appears in the radial 
contribution h—i to the backflow. It is this term which 
determines the single impurity contribution to the resid- 
ual resistivity ^,0. This is a simple example of the 
so called Das-Peierls theorem |],|l(J in electromigration, 
which states that the total force on the impurity is pro- 
portional to the additional resistivity it causes. The fact 
that the dipolar contribution ho to the backflow does not 
contribute to the friction force can be understood most 
easily by considering the linear response regime. Indeed 
to linear order in V the response at low velocities is purely 
reactive ||, while a finite resistivity can only appear at 
order V 2 . More generally, the coefficient ho is odd in 8i, 
while the force must be an even function of the phase 
shifts. This situation is in fact very similar to the case 
of a classical, incompressible and viscous liquid. Calcu- 
lating the backflow current around a sphere of radius R 
with boundary condition v = at the surface, one finds 
that 5j{x) has a contribution proportional to 1/r and 
a dipolar one. The associated function h(q) as defined in 
(1) is thus of the form 

h_ 2 

lim h c i(q) = —5- + h + . . . . (17) 

The coefficient of the 1/r contribution is h-2 = 6nRn 
while the strength of the dipolar backflow is negative and 
given by h Q — —irR 3 n (the corresponding problem in two 
dimensions has no solution which is known as the Stokes 
paradox). Calculating the associated friction F = —r/gv 
in a fluid with kinematic viscosity v it turns out |jj that 
only the leading term /i„2 contributes to r/s = 6nRn ■ rav 
while the dipolar backflow again drops out. Comparing 
the Stokes result with that for a Fermi liquid, we see that 
the fermionic friction coefficient for a scattering potential 
with characteristic range R such that a tr = irR 2 is equal 
to that of a classical liquid with finite kinematic viscosity 
vp = VpR/6. With typical values R = 2 A and vp — 
1.5 ■ 10 8 cm/sec for electrons in metals, we obtain vp = 
0.5 cm 2 /sec which is about fifty times the viscosity of 
water. From this point of view therefore, electrons in 
metals behave like a rather viscous liquid indeed ! 

In summary we have calculated the backflow around 
a fixed impurity in a Fermi liquid at low velocities and 



zero temperature. The dominant contribution is radial 
and decays like l/r d_1 . It is directly proportional to 
the current induced density modulation first discussed 
by Landauer and is responsible for the frictional force 
i.e. the finite resistivity. The subleading dipolar contri- 
bution of Pines and Nozieres has been evaluated to ar- 
bitrary orders in the scattering. It has been shown that 
the condition of perfect screening of a charged impurity 
gives rise to a novel sum rule for the corresponding phase 
shifts. We have also evaluated the total force on the im- 
purity. In agreement with the Das-Peierls theorem it is 
proportional to the scatterers contribution to the resis- 
tivity. For our simple situation this is just a consequence 
of Newtons third law. The fact that the dipolar term 
in the backflow gives no contribution to the force shows 
that the so called direct force in electromigration theory 
||,[ll| vanishes in the present case. This is a consequence 
of the way we have set up the problem: Instead of cal- 
culating a current as the response to an external field we 
have specified the incoming current which gives rise to 
a certain backflow or potential distribution It is an 
interesting future problem to develop a theory of electro- 
migration from this point of view, including the lattice, 
background scattering etc. 
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